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Abstract We study the geometrical influence on the Brownian motion over curved manifolds. 
We focus on the following intriguing question: what observables are appropriated to measure 
Brownian motion in curved manifolds? In particular, for those d-dimensional manifolds em- 
bedded in W^^^ we define three quantities for the displacement's notion, namely, the geodesic 
displacement, s, the Euclidean displacement, 5R, and the projected Euclidean displacement 
5R^. In addition, we exploit the Weingarten-Gauss equations in order to calculate the mean- 
square Euclidean displacement's in the short-time regime. Besides, it is possible to prove exact 
formulas for these expectation values, at all times, in spheres and minimal hypersurfaces. In 
the latter case, Brownian motion corresponds to the typical diffusion in flat geometries, albeit 
minimal hypersurfaces are not intrinsically flat. Finally, the two-dimensional case is empha- 
sized since its relation to the lateral diffusion in biological membranes. 
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1 Introduction 

Brownian motion occurs as a representation of a plenty of phenomena arising in various con- 
texts ranging from particle physics [IJ, general relativity [2] and condensed matter |3 1. In the 
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last decade there has been much interest in the study of diffusive processes on manifolds, mo- 
tivated by problems coming form biophysics |4 |. The transport processes ocurring on a biolog- 
ical cell is an interesting and complex problem. In particular, the motion of an integral protein 
through the plasma membrane has been approached from different point of views, which are 
themselves complementary. The most basic viewpoint is based on the Smoluchowski's equa- 
tion of a punctual particle on the membrane considered as two-dimensional regular curved 
surface. Further approaches consider membrane's thermal fluctuations ||ll-|l6l, as well as dy- 
namical fluctuations coupling to the stochastic motion of the protein f7^, and finite-size effects 
of the protein [8 |. 

In this paper we use the Smoluchowski's approach to study the geometrical component 
on the diffusion processes of the integral protein on the membrane. In general, we focus in 
the intriguing question: what observables are appropriate to measure Brownian motion in the 
present context? By analogy with Euclidean spaces it is not difficult to realize that the actual 
motion of the particle is through geodesic displacements, thus, the geodesic distance, s, is the 
proper notion for displacement on curved spaces ID, however, this quantity would be difficult 
to measure in an experiment. Nevertheless, unlikely Euclidean geometries for those curved 
manifolds embedded mR''+^ we have various quantities that also undergo stochastic dynamics, 
for instance, the Euclidean displacement 5R, and the projected Euclidean displacement 5R^. 
In this way, albeit Euclidean displacements are not the physical displacements they can be used 
as observables for the Brownian motion. Thus, from a practical point of view it would be easier 
to measure the Euclidean displacements than the geodesic displacement at least for Brownian 
motion on membranes (perhaps for the Brownian motion on curved space-time it is in the other 
way around because we are indeed immersed in that manifold). 

The Euclidean displacement has been, already, considered in ifTOl for Guassian polymer 
wrapping curved interfaces and also in llll to approach the transport modes on a cell mem- 
brane. The projected displacement has also been used in (5] and f6l. It is interesting that 5R^, 
as well as the Euclidean displacement, is refering to the ambient space where the manifold is 
embedded. In this sense, the mean-square (MS) geodesic displacement, s, will have an influ- 
ence from the intrinsic geometry whereas the Euclidean displacements, 5R and 5R±, will have 
an influence from the extrinsic geometry. In particular, in this work we study the mean-square 
(MS) values for the Euclidean and projected displacements using the method developed at 1121 . 
Also, this method allow us to give closed expressions for these mean values for spherical and 
minimal hypersurfaces for all time values. It is shown that the diffusion on minimal hypersur- 
faces measured from the observable 5R corresponds to the typical diffusion in agreement with 
the cubic minimal surfaces already studied in [JJJ and |[T3]| . 

This paper is organized as follows. In section 2, we summarized geometrical concepts used 
to describe intrinsic and extrinsic observables. In section 3, we present the diffusion equation 
and the displacement observables on curved marufolds and submanifolds. In section 4, we 
present the operator method used to evaluate the expectation values of observables. In partic- 
ular, we study the short-time regime of the observables 5R^, 5R and 5R^. In addition, we 
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show exact results for (5R^) and (5R) for spheres and minimal hypersurfaces valid at all 



time values. Finally, in section 5, we summarized our main results and we give our concluding 
perspectives. 

2 Preliminaries and notation 

In this section we review the preliminary notions about manifolds and sub-manifolds (following 
1141 and 1151 ') needed to describe the observables for Brownian motion on curved manifolds. 
Let M be a t/-dimensional manifold and C/ C M a local neighborhood. By definition of manifold 
U is locally diffeomorphic to a piece of Euclidean space. In particular, we are interested in 
manifolds endowed with a Riemannian metric g : Tp (M) x Tp (M) — M given by g = g„i, dyf ® 
dx'^, where gab is the meric tensor and Tp (M) is the tangent space for each p eM. Also, here, 
the Riemann tensor is denoted by R = R^b^^ ea®dxf' ®dx^ ®dxf^ and the Ricci scalar curvaturej^ 
by Rg. In addition, it is convenient to introduce the Laplace-Beltrami operator on scalars defined 
by Ag : & (M) -i- M given by 



where g = det g„i, and g is the inverse metric tensor. Also, the derivations are defined by 
da = d /dx", where xf with a = 1 , • • • , are local coordinates of some patch in the manifold. 

Let be a li-dimensional sub-manifold orientable in M''+^ with embedding functions 
X : C K'' — ;> Z C IR''+' which assigns (mi , • • • , m^) — > X (mi , • • • , m^). Each vector of the tan- 
gent space Tp (Z) can be spanned by {e^}, where ea := ^^X are the tangent vectors. The 1" 
Fundamental Form of the submanifold is defined by 1 : 7)^ (Z) — > M given by I (v) = v • v; thus 
the metric tensor adopts the simple structure gab = ^a'^b- Here, • is the canonical inner product 
ofK'^+i and I • I is the norm comming from this inner product. 

The normal direction to the tangent space is determined by the Gauss map N : Z — > 5'' 
defined by = 1 and N • eo = for each a. The curvature of the hypersurface can be under- 
stand in terms of the change of the Gauss map, thus the 2"'' Fundamental Form is defined by 
II : Tp{E) -^R given by II (v) = dNp ■ v; here the components of this form are the extrinsic 
curvature tensor Kab = &a • ^qN. The trace of this tensor is the mean curvature K = g"^Kab- 
Also, it should be noted that the "egregium" Gauss theorem implies that the Riemann tensor, 
Rabcd = KacKbd — KadKbc, depends just on the intrinsic geometry. The tangent space changes 
direction for each point p in the hypersurface. The manner how this change happens is captured 
by the Weingarten-Gauss structure equations 



The set {e,,} is a basis for the tangent space and is the corresponding basis in the dual tangent space. 




(1) 



V^efe = -Kab^, 

VfoN = K^"^a, 



(2) 
(3) 
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where is the covariant derivative compatible with the metric gab- 



3 Diffusion equation and observables on curved manifolds 

In this section, we introduce the simplest model to study Brownian motion on curved man- 
ifolds. This is a direct generalization of the diffusion equation on Euclidean spaces, which 
basically consist of replacing the Euclidean Laplacian by the Laplace-Beltrami operator Ag. 
This operator is often used to describe how a substance diffuses over a curved manifold. Also 
we may think it as heat diffusing on manifolds or polymer confined to curved surfaces IIIOI . 
Furthermore, it can be studied to determine the quantum propagator of a free particle on curved 
spaces 1 16 |. 

For a single particle diffusion over a manifold, we are interested in the probability density 
P : M X M X M+ ^ K+ such that P{x,x',t)dv means the probability to find a diffusing particle 
in the volume element dv when particle started to move at x' and it has passed a time t. This 
distribution is the same as the so called heat kernel as well as the propagator for a quantum 
particle moving on the manifold. In particular, here, the diffusion is governed by the diffusion 
equation (for P) 



dP{x,x' ,t) 
Jt 



= DAgP{x,x',t), (4) 



where D is the diffusion coefficient. Also, we required that P satisfies the initial condition at 
time t -5- 

\imPix,x\t) = —5'' ix-x) . (5) 

In addition, the existence of P (on a complete Riemannian manifold) is guaranteed if it also 
satisfies 

P{x,x',t) = P{x',x,t) (6) 

P(x,x',t) = [ dVyP{x,v,t~s)P(y,x',s) (7) 
Jm 

for any s > 0. This is the content of the theorem by Schoen and Yau |17|. These properties 
are easier to understand using the quantum analogy. The symmetry property (j6]l means that 
the evolution of a quantum particle from x to x' is the same evolution from x' to x, and the 
convolution property ^ means that this propagation from x to x' can be decomposed into 
individual propagations from x to y and y to x', for each value of y G M. 

In mathematics one of the challenges consist of estimate expressions for the heat kernel or 
the probability density (see II18I for a review) and one reason is because heat kernel encoded a 
strong connection with topological invariants of the manifold lf20l . For most geometries, there 
is not a closed form of the probability distribution. However, for short times there is a formal 
series solution for P(x,x',f) in terms of the Minakshisundaram-Pleijel coefficients [21], which 
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depends on both x and x' [.221 ; this series expansion is also called the parametrix expansion ||23l . 
It is noteworthy to mention that short-time depends on the especific geometrical dimensions of 
the manifold; below we will give a precise definition of what we mean by short-time in the 
present context. In what follows, we are interested in physical observables in order to have 
information about Brownian motion on manifolds. 

In order to get some insight about Brownian motion we often look at mean values of dis- 
placement and square displacement. These physical observables give us information of how 
particle diffuses in the space. For the Brownian motion on curved manifolds we have, in addi- 
tion, other random variables that are useful to understand the phenomena. These quantities are 
continuos scalar functions : M — > M defined on the manifold and, normally, are functions in 
C°° (M), although this is not the general case. The only rigid condition for an observable (x) 
is that its expectation value respect to P, defined in the standard fashion 

{i?{x))= I dvi?{x)P{x,x',t), (8) 
Jm 

is well-defined for all points x' in the manifold and for all time values. Note that {^{x}) depends 
on the initial position x'. Since we are concerning about physical observables for Brownian 
motion on curved spaces, we would like to address the question about what functions (jc) on 
the manifold are useful to describe Brownian motion. In particular, we refer to the displacement 
of a Brownian particle on the manifold. 

Let us recall that a free Brownian particle over an Euclidean space undergoes statistical 
fluctuations that are isotropic and homogenous. These properties imply rotational and transla- 
tional symmetries of the Euclidean space. In addition, these symmetries appear already in the 
Laplacian of the diffusion equation. Furthermore, the probability density inherited this invari- 
ance. For this case, the displacement is one observable of interest. The displacement is given by 
5R = X — X', where X and X' are two vectors in the Euclidean space; X' is the started point of 
the motion of the particle. Using the probability density for this case it is not difficult to show 
that (5R) = 0, and (5R^) = 2dDt; indeed, these observables capture the rotational and trans- 
lational invariance. Now, the Brownian motion of a free particle over curved manifolds will 
inherit the symmetries of the manifold. The displacement in this case is given by the geodesic 
distance, defined as follows. 

Geodesic displacement. In the intrinsic point of view, one particle displaces from p\ io p2 in 
M throughout a differentiable curve y : / C M — ^ M. Thus the geodesic displacement (GD) 
or geodesic distance, i^i = i, is defined to be the infimum length of geodesic curve beteween 
y(fi) =Xoj3(fi) andy(f2) = Xo j3 (fj), that is 

s = mi r\Y'{X{t))\dt, (9) 
Jti 

where X : ^ c M'^ — > M is a parametrization of the manifold and j3 : / C M — > M'^ is the pre- 
image of curve y. The expectation value, (i^), captures geometrical data of the manifold, and 
it gives the manner of how the intrinsic geometry influences Brownian motion. 
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In addition, for those li-dimensional sub-manifolds, embedded in M''+\ we have Euclidean 
displacements defined as follows. 

Euclidean displacement. In the extrinsic point of view, for Brownian motion over d-dimensional 
hypersurfaces (C there is another observable that is referred to the ambient space M''"'"^. 

The Euclidean displacement (ED) is defined by 5R : M — > IR''+^ given by 5R = X — Xq for 
all points X G M, where X is itself a parametrization of the hypersurface. This quantity also 
describes the displacement of the particle, but seeing it from the ambient space. Let us denote 
^2 = and i^?3 = |5Rp, where the distance function is |5R|. 

Projected displacement. Brownian motion also shows peculiar features in sub-spaces of the 
ambient space. In particular, here, we are going to consider only projections of the hyper- 
surface. The projected subspace is defined by 7r(M''+') = M'^, where the projection map, 
K : ]R''+' — > K'', is defined as usual by 7r(V) = v for V = (v,vo) e ]R''+'. In this subspace, 
the PROJECTED DISPLACEMENT (PD) is defined by the composition map 5Rx = no 5R. In par- 
ticular, for a domain of M such that it can be covered with one coordinate neighborhood we 
are able to use the Monge parametrization X = (x,/j(x)), where h{x) is the height function 
for X G W'. In these terms we have 5R^ = x — xq. Let us denote ^4 = (5R^)^. Note that this 
observable is also an extrinsic measure of the displacement of the particle. 



4 Expectation values of observables 

The general problem is to find the mean values for ^1, ^2, ^3 and ^4 for an arbitrary manifold 
and sub-manifold. In principle, these expectation values can be evaluated through the formal 
series solution of the diffusion equation in terms of the Minakshisundaram-Pleijel coefficients 
mentioned above. Here, we use an operator method introduced at |[T2ll . This method is inspired 
in the original calculations made by Perrin in his seminal papers about Brownian motion on 
spheres |24|. Next, P is a density distribution satisfying (|4|, Q, (|6]l and (|7]l. The method is 
encoded in the following result. 

Proposition 1 Let ^ -.M^R be an observable either in C°° (M) or W'+'^ x C°° (M), then the 
expectation value of 0, with respect to a probability density P, have the following derivations 
respect to time 



D'' j dv A^g0{x)P{x,x',t)+D^ J dvS/aJ", (10) 



where 



(11) 
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The proof of this proposition is along the following lines. First, let us differentiate {^{x)) 
with respect time, then substitute the diffusion equation Q. Next, we use the Green formula 
(45 1 and substitute the initial condition (|5]l. We gave a proof in 1121 . but this same result is also 
proved in II18I . 



Remark 1 Note that for each kis a. vector field on 



thus, by the divergence theorem, for 
compact manifolds the right hand side of Eq. ( 10 1 vanishes f251, except for non-trivial topolo- 
gies like circle or torus T^; these cases will be analysed elsewhere in [27] . For manifolds 
with boundaries we choose that P and VP vanish at the boundary therefore right hand side of 



Eq. ( 10 1 also vanishes. In particular, for these cases we have that 



(=0 



(12) 



In what follows, we will consider compact manifolds or manifolds where P and VP vanish 
at the boundary. In addition, assuming that {^{x)) /^f*^|,_Q are well defined on for a 
given M, we define the remainder R„ (t) by 



(^W) = I^W+/?«(0, 



(13) 



where the terms Gf = A^i^j are purely geometric factors. Thus by Taylor theorem (jsjl 
the remainder can be written in terms of the expectation value itself. The definition of the 
remainder, equation ( 13 i, is useful if we are able to prove that lim„_j.oo/?n (0 = because in this 



case 1^6 (x)) has a series Taylor representation 



k\ 



{Dtf 



(14) 



It is notable that not all expectation values satisfies lim„_j.„o/?„ {t) = 0; in [27 1 will be an exam- 
ple of this. For a given observable the difficulty lies to evaluate the terms Gf and to estimate 
lim„_i.<x,/?„ (f). The equation ( 13 I is very useful to access the short-time regime of the Brownian 
motion for the general manifold case, but also it can be used to find closed formulas valid for 
all times for some specific manifolds. 



4.1 The Brownian motion at short-time regime 



In what follows, we are going to give some estimations for the mean values of the observables^i, 
i^2, ^3 and ^4 at the short-time regime, that is for times t ^ Tg, where TG = 3d/ \Rg | D is called 
geometrical time |fT9]| . Basically, we use the proposition (jlj) and we compute the factors G^ for 
these observables with k= 1,2,3. 
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4.1.1 Intrinsic observables on manifolds 

The mean-value of s^ will capture intrinsic geometrical data of the manifold and it will give how 
this geometry causes a change in the standard diffusive behaviour. The geometric factors G^' 
cannot be written, in general, in a closed form for each k. However, in 1121 we have shown a 
formula for the mean-square geodesic displacement for the first three values k= 1,2,3. Hence, 
this value can be written as 



(s^)Ki2dDt- 



R,{Dt) 



1 

3! 



15 



nabn 



16 



abed 

45" 



nciDca fn 
K \Kdbca 



R. 



dcba ) 



1 



8"%] Rat 



(Dtf 



(15) 



where w is defined through the theorem Q. This result shows how the mean-square GD is 
deviated from the planar expression by terms which are invariant under general coordinate 
transformations. As a consequence of the Gauss "egregium" theorem all these terms are iso- 
metric and are built with 0(d) invariant combinations of the Riemann tensor. In principle, this 
result is valid for every Riemannian manifold endowed with metric tensor gai,- Clearly, in a lo- 
cal neighborhood Brownian motion is not affected by the geometry of the manifold, but as far 
away as the particle reaches the boundary of this neighborhood the curvature effects become 
apparent. 

Example 1 On spheres S'' , where Riemann curvature is Rahcd = ji {gacghd — gadgbc), the ex- 
pectation value of s^ is given by 



2dDt - ^ {Dtf + - ^ '- {Dtf + ■ 



2d{d~\) 
3 /?2 



45 /?4 

This result has also been obtained by direct calculation at |[T2ll . 



(16) 



Example 2 For regular surface embedded in 



Riemann tensor components are Rated — 
-2 {gacgbd ~ gadgbc), therefore the expectation value for these surfaces at short-time regime 
is given by 



4Dt 



15 



-K^G + 2AgKG 



(17) 



where Kq = 2Rg is the Gaussian curvature. Note that for developable surfaces. Kg = 0, the 
mean-square geodesic displacement behaves like the typical diffusion in this short-time regime 
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4.1.2 Extrinsic observables on submanifolds ofW'^^ 

For the Brownian motion over hypersurfaces of dimension d (submanifolds of we are 

interested to know the expectation value of 1^2 = I^Rp- In this case, also the geometric factors 
for ^2 cannot be evaluated in a closed form for a given integer k. Again, we are going to 
calculate the factors ^ for k= 1,2,3. For instance, for k = \ we have Gj ^ = '^g l^^j = 
2Vfl (5R • e''). Recalling the trace of metric tensor, g'^^ = d, and using the Weingarten-Gauss 
equation, ([3]), we get 

Gf2 = (2t/-2/f5R- N) 1 5R^o = 2t/ (18) 
In a similar way, by straightforward calculation, we obtain 

{\K)H ■ 5R - 2KS/" (Kj'eb) • 5R - 4 {S/^K) Kj'et ■ 5R - 2K^ 



G ^ - 



5R=0 



(19) 



and using geometrical identities in appendix ( |A.3[ ) we get 

Gf2 = 2K^KabK°^ - 2KAgK - 2Ag [K^) - 4Vfe (^KV"K/ + {V''K)Kf^ . 



Hence the mean-square Euclidean displacement (5R ) is written as follows 

{5^^)K2dDt-K'^{Dt f KW{K)+Ag{K^)+2VhJ''iK) {Dtf- 



(20) 



(21) 



where the scalar f [K) and vector J" (K) are defined as 

W{K)=AgK-KK,hK'"', 

J" {K) = KVhK'"' + 2K"^VbK. (22) 

This result shows how the mean-square ED is deviated from the typical diffusion behaviour. 
This deviation is also given by terms invariant under general coordinate transformations, but 
now they are referred to the ambient space where hypersurface is embedded; they are also in- 
variant under global rotations and traslations on the hypersurface. Like the geodesic displace- 
ment, in a local neighborhood, mean-square Euclidean displacement reproduces the standard 
Einstein kinematical relation. Also, when particle reaches the boundary of the neighborhood, 
the curvature effects emerge. For this case all curvature terms are built with -invariant of 
the second fundamental form or the extrinsic curvature tensor Kgi,. In other words, this observ- 
able encoded extrinsic information of the hypersurface. 

Both observables, and (5R)^, reproduce the standard mean-square displacement at very 
short-times, t ^Tq. It seems intuitive that at these short- times there are not at all curvature 
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influences on the Brownian motion because the local neighborhood looks like Euclidean space. 
However, even for those times the xobservable ^3 shows curvature effects on the Brownian 
motion. Following the same procedure for A; = 1 , 2, 3 we find 

(5R) w -KNDt - ^ N + /efc {Dtf 

(23) 

Therefore for non-zero times ? <C Tg = 3/2 | Z) there is still a contribution from the curvature 
on (5R). At short-times the normal direction of (5R) is explained as follows. Since there is 
not any preferential direction, tangent components of 5R cancel out in average and then, by 
symmetry, (5R) is along the normal direction. Nevertheless, as soon as the particle reaches the 
boundary of the local neighboorhood tangent components also contribute. 



Remark! For minimal hypersurfaces embedded in the mean curvature K is zero thus 
W (K) = and J" (K) = then the expectation values (5R) and (5R^) appear to be the same 
as those for the Brownian motion on flat spaces whereas expectation value of shows an 
influence of the curvature, since Rg < 0. 

Example 3 1 For a c?-dimensional sphere 5^ of radius R, the components of the second fun- 
damental form are Kab — ^gab^ thus the mean curvature is K — d/R. Then, we have ^{K) — 
—cP'/R^ and J'^iK) = 0. The started point in this case is Xq = RN. Therefore the expectation 
values of SR and 5R^ are 

(5R2) « 2dDt - ^ (Dtf + ^ {Dtf + ■■■. (25) 
This shows a different behavior in comparison to the one found using observable s. 



4.1.3 Extrinsic observable in the projected subspace 

In the subspace W'-, the projected displacement was defined by 5R_l = K (X) = x, for xq = 0. In 
order to compute the mean-square PD, (SR^^), let us write 1^4 = SR^^ = &2 — h^ and 5R_l = 
X — ftk, therefore the geometric factors satisfy 



Gf = Gf - A^ 



(26) 

x=o 
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Thus, for the geometric factors with k = 1,2,3 we have to calculate Api^ and A^h^- 

For k= I, we have Agh^ = IVahVh. Using the expression for the metric in this parametriza- 

tion, latter factor can be written as A^h^ = 2 ^"^'l! and using the normal vector N in this 
parametrization we find 



i2u2 „„A /I 3 7,3 



A„h^=2(l-N^). (27) 



The terms A^h and Apt are left expressed in covariant form. By straighforward calculation 
we find 

= 4 (v«V^/t) (V«Vfo/r) +4 (V)' + 12 (4gV„/i.) (V«/z) , (28) 

and 

4|/j2 = 8 (v«v''V^/i) (V<,VfoV,/i) + 32 (zigV''V''/j) ( V„ V,,/i) + 324^ (V«/j) (V"/j) 

+ 104g {Aghf - 8 (Ag/z) . (29) 

Note that V^, is the covariant derivative compatible with the metric gat and therefore it itself 
depends on the height function h. As in the previous cases, the mean-square projected displace- 
ment is written as 

(5Ri> ^2{d-l +N^) Dt + ^Cf (Dtf + ^Gf' (Dtf + ■■■ (30) 



For very short-times / <C Tg, mean-square projected displacement has the typical diffusion 



behavior in flat geometries (5R^ ) = 2dDprojt, but with a new diffusion coefficient 



Dproj = ^{d-l+N^), (31) 

modified by a determined geometrical content. This means that for the Brownian motion ob- 
served from the projected sub-space the diffusion is reduced, since Dproj is smaller than D. 
This modification is just a geometrical effect due to the point of view from where Brownian 
motion is seen. However, we can always rotate the hypersurface such that = I at x = Xq; 
after such a rotation we get Dproj = D. 



Remark 3 The result ( 3 1 1 has been obtained by different methods for two-dimensional surface 
at JSl and ||6l within the context of lateral diffusion of integral proteins in biomembranes. In 
these works they also consider the thermal fluctuations of the membranes. For instance at 0, 
under the basis of the Helfrich-Canham model ll26l for fluid membranes it is computed the 
effective value of the diffusion coefficient when thermal fluctuations are considered. It would 
be interesting to evaluate the next contributions of order |27|. 
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4.2 The Brownian motion for all time values 



We consider now the whole series ( 13 i for a particular set of observables. For these obser- 



vables, as we shall see, we are able to give exact and closed formulae valid for all time values. 
Let us start we the following 

Proposition 2 Let G (x) be an eigenfunction of Laplace-Beltrami operator Ag with eigenvalue 
—X, then the expectation value of G (x) is given by 

(^(jr)) = ^(jr')exp(-AZ)0 (32) 

Proof It is clear that 6 (x) fulfill all conditions: 6 (x) is a differentiable function. Indeed, the 
^-th action of is given A^^G(x) = {—X'fGix). The remainder R„ (f) have the following 
expression 



i-XD) 



/' 

Jo 



dT{G{x)){t-T)"dT, 



(33) 



then we have 



l^«(f)l < 



i-XD) 



n+i 



dx{G{x)) 



fit-T)" 

■Jo 



M 



-XD) 



n+l 



(34) 



(n+l)! ' 

where M is a number independent of n. It is elemantary that for any number a and £ > we 



have < e for sufficently large value of n, therefore lim„^oo/?„ (f) = 0. Now, using ( 13 i 

we get the wished result (|32]). 



□ 



Proposition 3 Let G {x) € C^^' (M) such that A^G{x) = C for each point on M, where C is a 
non-zero real constant, then the expectation value of G [x) is given by 

'x)) = G{x') +CDt 



(35) 

Proof It is clear that G{x) fulfill all conditions of proposition ([l]i: ^(x) is a differentiable 
function. Indeed, the k-th action of A^ is given by A^G{x) = for > 1. In this case, the 
remainder satisfies lim„^oo Rn {t) = 0. Now, using ( 13 i we get the wished result ( 35 i. □ 
Motivated by last two propositions let us open the following questions. Let X be a para- 
metrization for an Euclidean sub-manifold and let s be the geodesic distance in a manifold. 
Last result leads us to pose the following questions for X real, 1. What is the collection of d- 
dimensional sub-manifolds of such that each submanifold has at least one parametrization 
X satisfying each of the following conditions (i) AgX = XX or (ii) AgX^ = AX^ or (//;') A^X^ = 
A ? 2. What is the collection of t/-dimensional manifolds such that each of them satisfies the 
following conditions (/) AgS^ = A or (jj) A^s^ = Xs^ in at least one local neighborhood? Answer 
of first question with condition (;') is given by the following two propositions. 
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Proposition 4 The submanifold M, embedded in M'^+^ is a d-dimensional sphere S'^, with 
radius R, if and only if there is a parametrization X such that — ^^X = -^\,for non-zero R. 

Proof On one hand, let us assume that the submanifold M is part of a li-dimensional sphere S^, 
with radius R, then there is a parametrization X such that X = 7?N, where N is the normal vector 
point outward the hypersphere. Also, the mean curvature oiS'' isK = d /R. By the Weingarten- 
Gauss equations AgX = ~KN, therefore — 4gX = |^X. On the other hand, if —AgX = AX 
by the Weingarten-Gauss equations we get X = that is X^ = therefore V,, (K^ /X) = 
X • Cfl = ^ A' is constant <^ M is 5"'. □ 



Proposition 5 The submanifold M, embedded in Mf^^^, is a d-dimensional minimal hypersur- 
face, K = 0, if and only if there is a parametrization X such that ~Ag\ = 0. 

Proof On one hand, let us assume that the sub-manifold M is a c/-dimensional minimal hyper- 
surface, namely mean curvature vanishes, ^ = 0, then by Weingarten-Gauss equations there is 
one parametrization such that Ag\ = 0. On the other hand if there is a parametrization such 
that AgX = 0, then by the Weingarten-Gauss equations KN = 0, that h K = 0. □ 
We ignore the answer of first question for the conditions (ii) and (///). However, for condi- 
tion (//) it is clear that spheres are examples for X =0. Also, it is noteworthy to mention that 
condition (ii) is equivalent to the relation 2d — 2KX ■ N = AX^, thus minimal hypersurfacess 
do not belong to the collection of this condition. Also, minimal hypersurfaces belong to the 
collection associated to the condition (///) since for these sub-manifols we have AgX^ = 2d. In 
addition, we ignore the answer of question two, but for condition (/) it is clear that at least flat 
geometries are examples and for condition (jj) we do not know even if there is a manifold that 
satisfying AgS^ = Xs^. All, these results are useful to proof a general structure for expectation 
values of 5R and 8'S? for minimal hypersurfaces (K = 0) and spheres S"^. These results are 
encoded in the following theorems. 

Theorem 1 The expectation values of5K and 5R^, with respect to P, for each minimal hyper 
surfaces of dimension d are given by 

(5R) = 0, (36) 
(5R2) = 2dDt, (37) 

for all values of time t. 

Proof Since the sub-manifold is a minimal hypersurface then by proposition ^ we have 
Ag\ = and by proposition (j2jl we get (X) = Xq <^ (^R) = 0. For the mean-square ambient 
displacement, let X a parametrization of the minimal hypersurface then by the Weingarten- 
Gauss equations AgX^ = 2d, therefore for the proposition (jsj we get (X^) = Xq + 2c/Z)f, then 
(dR^)=2dDt. □ 



Remark 4 This theorem ([T]) is consistent with the general formula at the short- time regime (21 



see remark Q. This is in agreement with the result for the cubic minimal surfaces explicitly 
studied in fTTl and lfT3l. 
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Theorem 2 The expectation values ofdR and 5R^, with respect to P, for spheres S^' are given 
by 

(5R)=Xo(^exp(^-^Z)f^-l^, (38) 
(5R2)=27?2(^l-exp(^-^Z)f^^, (39) 

for all values of time. 

Proof Since the sub-manifold is a sphere 5"' by proposition ^ there is a parametrization X 
such that — Zi^X = ^X, that is X is an eigenfunction of A^, with eigenvalue — therefore by 

proposition (2i we get (X) = Xoexp (^—-^Dt^, where Xq is the started point. Now, since the 

sub-manifold is an sphere, thus AgX^ = 0, then (5R^) = (X^) + (X^) - 2 (X) • Xq, therefore 
by proposition (j2| and ([3]l we get the wished result. □ 



Note that equations ( [38| l and ( [39| l at the short-time regime reproduce ( [24| ) and ( [25[ ), re- 
spectively; see example Last result has been already found at [ 10| and [ 1 1 1 by alternative 
methods. Now, let us consider a d-dimensional infinite cylinder with radius R. This cylinder can 
be thought as Cyl = S'^~^ x K. In this case, the embedding functions can be written in terms 
of that {d — 1) -dimensional sphere X = (X^^-i where z G K. The metric can be written as 
gab = diag (Ijg/y), where gij being the metric of S'^~^ and the Laplace-Beltrami operator in this 
case is given by ^cyi = ^^rf-i + d^/dz^- The started point is chosen to be Xq = (1,0, • • • ,0), 
thus we have the following 

Corollary 1 The expectation values of 5R and 5R^, with respect to P, for cylinders Cyl are 
given by 

{5R) =Xo(cxp(-^Dt)-l], (40) 



(5R2) = 2Dt + 2R^ ^1 -exp ^-(n- 1) 



Dt 



(41) 



for all values of time. 



Remark 5 In particular, we have {5z) = and (5z^) = 2Dt which corresponds to the one 
dimensional diffusion. 

Example 4 Let us take a hemisphere of radius R and let i7 = be the projected subspace 
from this hemisphere. Thus a parametrization of this hemisphere is X : J7 — > defined by 



X(9,p) = ypcos(p,psin(p,y R^ - j . (42) 

This means that metric tensor comes in the following form gab = diag [p^ ,R^ / [R^ ~ p^)) and 
the projected displacement becomes dR± = (pcos(p,psin9,0). In particular, we can verify 
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that Ag (5R^) = 4 - for all values of p G M+, thus it is not difficult to show that by a 

straighforward calculation we find Al (5Ri) | = 4 f - ^ j for all naturals k^Q. Also, 

for this observable, 5R^, we can verify that the remainder 7?„ [t) goes to zero for large n. 
Therefore in this regime the mean-square projected displacement is given by 

(5Ri) = -/?2(^l-e (43) 
valid at all time values. For short-times this expression reproduces the general structure for 



3.5 



3.0 



2.5 



2.0 



Q 



1.5 



1.0 



0.5 





ADt I 
R- /' 


(^-4)/2: 






2 




,' '"^ 








2/3 









0.5 



1.0 



1.5 



2.0 



2.5 



Dt 



3.0 



3.5 



Fig. 1 Mean-square geodesic and Euclidean displacements as a function of time for a Brownian particle dif- 
fusing on the sphere. The expectation value (.s^) corresponds to our result in |27 1; the expectation value (5R^) 
corresponds to the equation |39| and the expectation value (5R^) corresponds to Eq. j43| . The straight lines 
stand for the short and long-time limits. 



the projected mean-square displacement (30 1. In particular, the diffusion coefficient does not 
change since A^- = 1 at p = 0. 

In figure we show a comparison between expectation values of {s^), (^R^) and (5R^) 
for all time values; these three observables coincide to ADt for short-times, however, for large- 
times it is shown that the mean-square PD saturates faster than mean-square GD and ED do 
as a consequence of the projection to the plane. It is also interesting that the saturation value 
(for large-times) goes to the value (5R^) = (2/3)/?^ rather than 2R^ as it happens for the 
observable /5R^V 
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5 Concluding perspective 

In this paper we have studied Brownian motion over curved manifolds and, mainly, on Eu- 
clidean sub-manifolds of dimension d. Our framework is based on the diffusion equation on 
curved manifolds. Here we have addressed the question about what functions G (x) are useful 
to described Brownian motion on curved spaces. In particular, we reviewed the notion about 
displacement and we define three quantities for this notion, namely, the geodesic displacement 
s, the Euclidean displacement, 5R, and the projected displacement, 5R^. It is noteworthy to 
mention that the controversy posed in ||9| about the displacement is figure out by consider that 
all these displacement capture information of the Brownian motion. Furthermore, from practi- 
cal point of view it would be easier to measure extrinsic displacements than the geodesic one 
because of the embedding of the membrane [6J. 

Here, we study the short-time regime of expectation values of 5R, 5R^, and 5R^ using the 
operator method introduced at lfT2ll . It is remarkable that for sufficiently small times, t <^ Xg, 
the quantities (5R^) and (5^) show no curvature effects while 5R shows them. Indeed, at 
this short-time regime the mean value for 5R is proportional to the scalar curvature K and 
it is pointing outward the hypersurface at the normal direction N. In other words, it is like if 
local curvature has been manifesting through a normal force acting on the particle. For the 
projected distance 5R^ it is found that diffusive coefficient D is modified by the projected one, 
Dpioj = ^ [d — \ +N}^, by purely geometrical effects. However, we point it out that for a fix 
geometry we can always find a global rotation such that D = D^^^y Nevertheless, this is not 
the case for membranes under thermal fluctuations because, in general, average of N'} is not 
going to be equal to one |6|. In addition, it is shown that expectation values of eigenfunctions 
of Laplace-Beltrami operator can be written in a closed form. Particularly, it is proved that 
expectation values are linear in time for functions that satisfy Agi^ = constant. These results 
allow us to find closed expressions for expectation values corresponding to spheres and min- 
imal hypersurfaces. In particular, we found surprising that for the latter case the mean-square 
Euclidean displacement is exactly the one found for flat geometries, although the Brownian 
motion has intrinsic curvature effects. A very similar effect happen for developable surfaces 
(with zero Gaussian curvature) at least for short-times [gl. 

For the future work, on one hand the surprise found for the minimal hypersurfaces using 
5R immediately generalizes to higher order momenta. Even, it opens the possibility of finding 
an exact result for the density probability where Weierstrass-Enneper representation may play a 
central role. On the other hand, the behaviour of the mean-square displacement (either intrinsic 
or extrinsic) for fluid membranes (i.e. described, for instance, by the Helfrich-Canham model) 
under thermal fluctuations can be study, relatively simple, at least to the one-loop order in 
Feynman diagrams. 
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A Some important theorems and identities 

A. 1 Boundary terms and Green formula 

The boundary terms in proposition ^ can be re-written as follows 

Remark that this expression involves the terms P and V"P, then by imposing the mixed Neuman and Dirichlet 
boundary conditions P|^m = ^^'^ ^"^\dM = we are able to ignore the second term of equation jlO^ . 
A key ingredient for the proof of the proposition (|T|l is the Green formula. Let (j>i and (j>2 two scalar function on 
the manifold M then the following identity is satisfied 

/ dv<piAg(l)2= I dv {Ag(l>i)(j>2+ I da(^idv<h-[dv^\)^), (45) 

Jm JM Jam 

where da is the volume element of the boundary, i?v = v • V and v is the outer normal vector of the boundary 

dm. 



A.2 One real variable theorems 

In this subsection let us follow |28|. Let / : / C K ^ K a function. Let us denote /<"' [t) = d"f{t) /dt" and 
ail = ■' ' . In addition, let us define the Taylor polynomia 

PnM{')="o + ai {t-a) + ---+a„{t-a)" . (46) 

Theorem 3 (Taylor theorem). Let us suppose that /',••• ,/"^' are defined and integrable in \a,t\ and that the 
remainder Rn,a {x) is defined by 

/(?)=/ («) +/' {a){t-a) + --- + (f - fl)" +/?,,„ it) , (47) 

nl 

then 

Theorem 4 Let / : / C K — > K a function for which f (a) , /"+' (a) with aEl exist, then 

^ m-P^(') ^0 (49) 
t^a [t — a) 

In the sense of this theorem for the limit where t ^ awe write / (f) ~ Pn.a (') ■ 
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A. 3 Geometry identities 

The following identities are useful for the calculations of the mean-values. This identities can be straighforward 
find them using the Weingarten-Gauss, |2j and ([3J, several times . 

Va im = {VaK) N + KKab e" (50) 
4yN = V,^:°*efc-rt„„N (51) 
Ag (^N) = {AgK - KKabK"") N + [kV.K"" + IK^'V^k) (52) 

VcAg(m) = |Vc (^AgK - KKabK'*'^ - (^KVaK"'' + 2K"''V„K^^>i 

+ I (AgK - KKabK"'''^ K/ + V,. {kV^K"'' + 2K"''VaK^ K"'' j ej (53) 

(54) 

In particular, identities ([50}-|53j are useful to determine the mean and mean-square Euclidean displacement 
5R. Now, for the projected displacement, 5Rx are useful the following identities. 

V„5R_L =e„-V„ftk (55) 
4j 5R_L = - (^N + Aghk) (56) 
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